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Dielectric tensor of monoclinic Ga 2 03 single crystals in the spectral range 0.5 — 8.5 eV 
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The dielectric tensor of / 3 -Ga 2 03 was determined by generalized spectroscopic cllipsometry in a wide spectral 
range from 0.5 eV to 8.5 eV as well as by calculation including quasiparticle bands and excitonic effects. The 
dielectric tensors obtained by both methods are in excellent agreement with each other and the observed 
transitions in the dielectric function are assigned to the corresponding valence bands. It is shown that the 
off-diagonal element of the dielectric tensor reaches values up to \s xz \ ss 0.30 and cannot be neglected. Even 
in the transparent spectral range where it is quite small (|£ X2 | < 0.02) it causes a rotation of the dielectric 
axes around the symmetry axis of up to 20°. 
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I. INTRODUCTION 


Transparent conductive oxides (TCO) are important materials for the realization of transparent electrodes 1 , e.g. 
for solar cells 4 and transparent transistors 5 . Indium tin oxide (ITO) which is up to now the material of choice in 
industrial applications has the drawback that indium is quite expensive and therefore much effort is devoted to develop 
indium free TCO 2 . One promising material is Ga 203 with a reported band band gap energy of about 4.8 eV 6,7 so 
that even in the presence of crystal imperfection high transmission in the visible and even in the UV-A/B spectral 
range can be obtained. Furthermore, the band gap energy is larger then the absorption edge of the atmosphere so 
that Ga 203 is also a promising candidate for solar blind photodetectors 8,9 . 

At ambient conditions, Ga 203 crystallizes in a monoclinic crystal structure (/3-phase) and the dielectric function 
is a tensor rather than a scalar function. Although the optical characterization of thin films was done by several 
groups 8,10-15 , rotation domains present in the thin films and the inferior crystal structure compared to a single 
crystal prohibited the determination of the tensor components and only an effective dielectric function was measured. 
Experiments on single crystals are rare and mainly limited to transmission 6,16 and absorption 6,11 measurements. 
From such data the entire dielectric tensor cannot be determined and only the absorption edges as well as the refractive 
index in the transparent spectral range were deduced. This strongly limits the comparison with theoretical calculations 
of the band structure. Besides experimental work, the band structure as well as the tensor of the dielectric function 
were calculated by several authors 7,18,19 . However, for the reported dielectric function the crystal was treated as an 
orthorhombic system, i.e. the dielectric tensor consists only of diagonal elements and the off-diagonal element was 
neglected. 

In this work we present the full dielectric tensor of P-G& 2 O 3 within the spectral range 0.5 — 8.5 eV. We show that the 
off-diagonal element {e X z) has a significant impact on the orientation of the dielectric axes. This knowledge is important 
for polarization sensitive applications. The paper is organized as follows: first we describe the experimental setup, 
the investigated sample and present the theoretical framework which we used for the determination of the dielectric 
tensor. In the following section, the obtained dielectric tensor and the rotation of the dielectric axes are discussed. 


II. EXPERIMENT AND THEORY 
1. Sample 

The equilibrium crystal structure of Ga 203 at ambient conditions is monoclinic, the so-called /3-phase. The unit 
cell of this structure consists of 8 Ga atoms which are surrounded by oxygen atoms where 4 Ga atoms exhibit an 
octahedral, the other four a tetrahedral coordination 20 . The C 2 h symmetry axis of this structure coincides with the 
crystallographic 6 -axis in real space and the lattice constants are about 1.22 nm, 0.30 nm and 0.58 nm for the a, b and 
c-axis, respectively. The angle 7 between the a and c-axis was determined to be 103.7° 20 . A schematic representation 
of the unit cell is shown in Fig. 1. 

For the investigations presented here, we have used two bulk single crystals (size of 10 x 15 mm 2 ) with different 
surface orientation, namely (201) and (010). They were fabricated by the Tamura Corporation by means of floating 
zone technique. The surface planes of the two crystals are marked in Fig. 1. The FWHM of the (102) X-ray diffraction 
peak was specified to be 89" (48") whereas the corresponding peak of the in-plane direction of the (102) plane is about 
46" (25") for the (010) ((201)) oriented sample. The as-received surface was investigated by atomic force microscopy 
(AFM) which revealed a smooth surface without atomic steps. The roughness was determined to be R s ss 0.2 nm 



FIG. 1. Schematic representation of the unit cell of P-G&zOz. The tetrahedrally coordinated Ga atoms are shown in blue 
whereas the octahedrally coordinated ones are shown in green. The oxygen atoms are marked in red. The (010) and the (102) 
planes are highlighted by the orange and red rectangle, respectively. (Image created by VESTA 21 ) 
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(0.1 nm). This smooth surface and the absence of the atomic steps indicates an amorphous surface layer with possibly 
different optical properties than the bulk crystal. 


2. Dielectric tensor 

The monoclinic structure of Ga 2 C >3 causes an optically biaxial anisotropic behaviour with a refractive index and 
absorption edge depending on the polarization. In the case of Ga 2 C >3 this dependence of the absorption edge on 
the polarization was explained by the symmetry of the valence bands and the conduction band by means of band 
structure calculation 19,22,23 . These calculations also yield, that most of the charges corresponding to the conduction 
band minimum are located in interstitial regions whereas in the case of the valence bands the charges are located at 
the p-orbitals of the O atoms. The width of the valence bands was estimated to be about 7 eV. Surprisingly, it was 
found that the conduction band minimum is nearly isotropic. Thus the anisotropy of the reduced mass is determined 
by the anisotropic behaviour of the valence bands 7 . 

The optical response of a material is described by its dielectric function (DF). For a monoclinic system this DF is 
a tensor given by 24 : 

( ^XX 0 X / \ 

0 £ yy 0 I • (1) 

^XZ 0 £ zz J 

The y-direction coincides with the symmetry axis. Since this is a symmetric tensor, it can be diagonalized for each 
photon energy separately and the corresponding diagonal elements are the semi principal axes of the ellipsoid of wave 
normals 24 . These axes of the ellipsoid are also often called dielectric axes. The dispersion of each element of the 
dielectric tensor is in general independent from the others so that the orientation of the dielectric axes with respect 
to the sample coordinate system depends on the photon energy 25 . This behaviour is called colour dispersion. For the 
tensor given by Eq. (1) one dielectric axis is fixed along the y-direction whereas the other two rotate as a function of 
the photon energy in the x — 0 -plane around the y-axis. 


3. Spectroscopic Ellipsometry 

Spectroscopic ellipsometry was applied for the determination of the dielectric function of 13 -G& 2 O 3 in the spectral 
range of 0.5 — 8.5 eV. This method determines the change of the polarization state of an incoming light after reflection 
on the sample surface 26,2 ' which is usually expressed by the ratio of the complex reflection coefficients for light 
polarized perpendicular (s) and parallel (p) to the plane of incidence. However, in the case of an optically anisotropic 
material, conversion of s- into p-polarized light takes place (and vice versa) and the change of the polarization state is 
then expressed by the 4x4 Mueller matrix M by 5 re f = MS- m with S- ln (5 re f) being the Stokes vector of the incident 
(reflected) beam 26 . For the interpretation of the data, the Mueller matrix is often described as a matrix of four 2x2 
block matrices. The diagonal block matrices are then mainly related to the different reflection coefficients for the p- 
and s-polarized light whereas the off-diagonal block matrix describes the mentioned conversion of the polarization 
from s to p (and vice versa). 

For the determination of the 4 independent complex quantities of the dielectric tensor (cf. Eq. (1)) the measurements 
have to be performed at different orientations of the dielectric axes with respect to the plane of incidence 28 . This was 
realized by using two bulk single crystals with the surface orientation mentioned above and by rotating each crystal 
around its surface normal by 0°, 45°, 90°, 105°, 180°, 225°, 270° and 305° in order to ensure a high sensitivity to each 
component of the dielectric tensor. The corresponding orientation of the sample system with respect to the laboratory 
system is then described by using the Euler angles (ip, i9, if>) in the yzx- notation 29 . For (ip, 1 }, -0) = (0,0, 0), the sample 
as well as the laboratory system coincide with each other. Since in Ga 203 the angle between the crystallographic a 
and c axis is non-orthogonal, there is no direction which defines a cartesian coordinate system in this plane. Therefore 
we choose in the following the system defined by: e x || a-axis, e v || 6 -axis and e z = e x x e y , e» being the unit vector 
in the i direction. As angles of incidence we have chosen 60° and 70°. 

For the determination of the dielectric function of Ga 203 from the recorded spectra we used a layer stack model, 
consisting of the semi-infinite substrate and a surface layer. Since both single crystals were fabricated by the same 
technique and produced by the same company, the dielectric function was assumed to be the same and analyzed 
simultaneously. The surface roughness of the samples was described by an effective medium approximation (EMA) 30 . 
For both substrates an effective thickness of about 2 nm was determined reflecting the smooth surface obtained by 
AFM. For the calculation of the DF of this layer, for simplicity only the diagonal elements of the dielectric tensor of 
Ga 203 were considered and mixed with the DF of void with a fraction of 50% : 50%. 
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A 

B 

(10“ 2 pm 2 ) 

C 

(10“ 3 pm 4 ) 

£xx 

3.479 

4.727 

1.753 

£ yy 

3.626 

4.739 

1.242 

£ ZZ 

3.574 

4.882 

2.129 

&XZ 

0.016 

-0.185 

-0.146 


TABLE I. Parameter of the tensor elements of the dielectric function of / 3 -Ga 203 in the transparent spectral range (E < 4eV). 


By applying this layer model, the DF was obtained by using a Kramers-Kronig consistent wavelength-by-wavelength 
analysis for each energy. By doing so, the DF exhibits a non-vanishing imaginary part in the order of 10 -3 in the 
visible spectral range. This causes a significant absorption in this spectral range which is physically wrong. This is 
probably caused by the unknown optical properties of the surface layer which can lead to such an effect, similar to 
an isotropic material where no surface layer is taken into account, and to the limitation of the used EMA model. 
However, due to the small thickness of this surface layer (4ma < 2nm) a determination of its DF, especially its 
anisotropic character, is quite challenging due to the low sensitivity to such a thin layer and the main results are 
unaffected by the choice of the surface layer. To overcome this problem, the spectral range was divided into two 
spectral ranges: a ’’transparent range 11 (0.5 — 4.5eV) where the absorption is set to zero and the ’’absorption range 11 
(3.0 — 8.5 eV) which includes the spectral range where absorption takes place. 

In the ’’transparent range 11 (0.5 — 4.0eV) the absorption was set to zero and the dielectric function can be fully 
described by its real part (eipj). For materials with an orthorhombic (or higher) symmetry in this case the relation 
rq = \[£~i between the components of the refractive index (n*) and the dielectric function holds where rii can be 
described by an Cauchy approximation. In our case, due to the presence of the off-diagonal clement this relationship 
is not valid anymore and we used 


£ i ,ij 


B 


v 

A 2 


C, 


_u 

A 4 ' 


( 2 ) 


Note, by using this approximation a negative dispersion can be achieved as it can be occur for the off-diagonal element 
of the dielectric tensor which is not possible by using an Cauchy approximation. The corresponding parameters are 
given in Tab. I. 

The second range includes the ’’absorption range 11 (3.00 — 8.65eV). Here the DF was obtained by applying the 
wavelength-by-wavelength analysis as described above. In the spectral range where both parts overlap (3 — 4eV) a 
good agreement between the DF deduced by the Cauchy approximation and the wavelength-by-wavelength fit was 
obtained. In order to obtain a single DF for the entire spectral range, the two parts of the DF were merged where a 
linear weighting function was applied in the overlap range. This weighting ensures a smooth onset of the imaginary 
part of the DF and was done in such a way, that the transparent (absorption) spectral range contributes by 100% 
(0%) to the overall DF at an energy of E = 3eV whereas at E = 4eV the contribution was set to 0% (100%). The 
overall dielectric function was than smoothed by using the Savitzky-Golay algorithm. 

We note that M. Dressel et al. proposed a bond polarization model 31 in order to analyse triclinic pentacene. In 
this model it is assumed that the dipole element and therewith the polarization is oriented along the crystallographic 
axes. The advantage of this model is, that the off-diagonal elements of the dielectric tensor can be expressed by the 
principle elements and the known angles between the crystal axes, which reduces the number of unknown independent 
tensor elements to three. Therefore, the major axes of the indicatrix (i.e. tensor related to the real part of e) coincide 
with those of the conductivity tensor (i.e. related to the imaginary part of e). Although this model was successfully 
applied to pentacene 31 and slanted columnar thin films 32 it cannot be applied in our case and in general since the 
major axes of the indicatrix and the conductivity tensor do not coincide with each other (see inset Fig. 3a). 


4 ■ Theory 

In a first step we determine the ground-state atomic structure by minimization of the energy with respect to all 
lattice parameters and ionic positions, employing density functional theory (DFT) 33 and the projector-augmented 
wave (PAW) method 34,35 as implemeted in the Vienna Ab-initio Simulation Package (VASP) 36,3 '. In order to obtain 
structural data very close to experiment we used the superior AM05 gradient-corrected exchange-correlation (XC) 
functional 38,39 . A plane-wave cutoff of 820 eV and a grid ofllxllxllk points has been used for the structural op¬ 
timization. The high plane-wave cutoff was necessary because we employ the calculated stress tensor for optimization 
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of the lattice parameters and this stress tensor may suffer from spurious Pulay stresses due to the incompleteness of 
the basis set 40 , demanding for high plane-wave cutoffs. 

The electronic excitations and optical spectra of /3-Ga203 are calculated employing state-of-the-art techniques 
based on many-body perturbation theory (MBPT) 41,42 . We use the GW scheme proposed by Hedin 43 for the XC 
self-energy, employing an HSE 44 hybrid functional as starting point for the electronic quasi-particle band structure 
to calculate single-particle excitation properties. The HSE starting point is chosen because the hybrid functional 
mimics already in a very crude manner parts of the quasiparticle effects. Hence the GW corrections remain much 
smaller than for a DFT starting point justifying much better a single-shot GoWo approach. For optical properties the 
inclusion of excitonic effects, treated within the Bethe-Salpeter equation (BSE) 41,42 framework, is crucial in order to 
allow a reasonable comparison with experimental spectra. We use the standard implementation of GW and BSE of 
VASP which employs the full frequency-dependent dielectric function without any approximations like the plasmon 
pole approximation 45 or model dielectric functions 46 for the calculation of the screened Coulomb interaction W and 
standard static screening in the BSE case 41 . 

For both, the GW and BSE calculations, we have to limit the k-point sampling to a grid of 8 x 8 x 8 k points 
including the Brillouin zone center T point. As a plane-wave cutoff we use a value of 410 eV. For the representation 
of the full microscopic inverse dielectric function entering the screened Coulomb interaction W in GW and BSE we 
can reduce the cutoff to 150 eV. All calculated optical spectra are smoothed with a broadening of 0.2 eV. Due to high 
computational demands we have to use a numerically very efficient time evolution scheme 4 ' to solve the Bethe-Salpeter 
equation where the calculation of the w-dependent polarizability can be considered as an initial-value problem 4 '. This 
avoids the diagonalization of a giant pair Hamiltonian matrix with typical ranks of 300.000 to 400.000. However, 
one looses information about eigenvalues and eigenvectors of the BSE pair Hamiltonian matrix. Anyway, if bound 
exciton peaks can be observed below the fundamental gap energy one can at least estimate the bound exciton binding 
energies from the peak positions relative to the fundamental single-quasiparticle gap. 


III. RESULTS 

5. Dielectric Function 

For selected crystal orientations with respect to the angle of incidence, the recorded as well as the calculated 
spectra of the Mueller matrix elements are shown in Fig. 2. A good agreement between the experimental and the 
calculated spectra can be recognized. As can be seen, the off-diagonal elements of the Mueller Matrix are strongly 
pronounced which reflects the anisotropic character of the sample and therewith mode conversion. This effect is 
especially pronounced in the absorption range. Only for the case, where the symmetry axis of the crystal (6-axis) 
is perpendicular to the plane of incidence, the off-diagonal elements vanishes and no mode conversion takes place 
(Fig. 2c). 

The DF derived from the best match between the experimental and calculated spectra is shown in Fig. 3a. In the 
visible spectral range, the largest anisotropy is observed between the tensor elements for the x and y direction to be 
E yy — e xx ps 0.15, whereas the difference between E yy and e zz is quite small (ps 0.05). In the absorption range, the 
situation changes. In the diagonal elements, several peaks are observable caused by the transitions from valence bands 
into the conduction band. The first transitions occur for the dipole polarized perpendicular to the symmetry axis, 
namely the z- and x-direction, at around E ps 4.8 eV. The first one for a polarization along the y-direction is shifted 
around 500 meV to higher energies. This causes a more pronounced dispersion in the transparent spectral range for 
the e xx and e zz component as for the £ yy . At two energies, namely E ps 3.26 eV and E ps 4.5 eV, the dispersion of the 
elements e zz and e xx crosses with the one of £ yy , respectively, causing two uniaxial points. 

One important property, which is as mentioned often neglected in the literature, is the monoclinic nature of /?- 
Ga203, i.e. the presence of the off-diagonal tensor element of the DF. In the transparent spectral range (E < 4eV) 
this element is small (|ei, X z| < 0.02) while in the absorption range several peaks are observable which reach values 
between —0.25 at E ps 4.8 eV and E ps 6.9 eV and 0.30 at £ ps 5.1 eV. As already mentioned, the presence of the 
off-diagonal element causes rotation of the two dielectric axes in the x — z-plane around the y -axis by an angle (f> in 
dependence on the energy and the optical property (i = 1, 2=real,imaginary part) given by 


tan 2^= 2£i ^ z .. (3) 

&i,xx ^i,zz 

Although the off-diagonal element in the transparent spectral range seems to be negligibly small, in conjunction with 
the dispersion of the diagonal elements, it causes a rotation by roughly 20° (inset Fig. 3a). 
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FIG. 2. Experimental (symbols) and model calculated (lines) spectra of the Mueller matrix elements for an angle of incidence 
of 70° for different orientations of the crystal (indicated by the Euler angles (tp, 9, i/j) on top of each column). 
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FIG. 3. (a) Experimentally determined (left column) and theoretically derived (right column) tensor components of the 

dielectric function, inset: Rotation angle of the dielectric axes around the y-axis for real (black) and imaginary (red) part, (b) 
The calculated band structure of /3- Ga 2 03 . The colour of the valence bands represents the polarization of the dipole-allowed 
transition into the lowest conduction band at the F-point (marked by arrows in (a)). The valence bands which have weakly 
dipole-allowed transitions are shown in light grey. 
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6. Band structure and theoretical dielectric function 


The band structure calculated for low temperatures (T = 0 K) in the vicinity of the upper valence bands and the 
lowest conduction bands is shown in Fig. 3b. By taking into account the symmetry of the bands and the corresponding 
transition matrix elements, the dielectric function is obtained and shown in Fig. 3a, being in excellent agreement 
with the experimental one. The additional peaks in the theoretical derived DF compared to the experimentally 
determined can be attributed to artefacts due to the limited k-point sampling used for the calculation and the small 
broadening of the transitions considered in the calculations. Beside the additional peak structure, an energy shift of 
about 500 meV between the calculated and experimentally observed transition is observable. Since the experimental 
observed transitions are blue shifted compared to the theoretical calculated, this shift can hardly be explained by 
the different temperature used for the calculation and in the experiment. A reason might be that electron-phonon 
interaction has to be taken into account in the calculations or maybe also additional screening effects due to free 
carriers, although a characterization of the sample indicates free-carrier concentrations below 10 19 cm -3 . Both might 
reduce the excitonic effects yielding a blue shift of the calculated spectra. Also strain in conjunction with a rather 
large deformation potential can be a potential reason. However, our calculated lattice parameters a = 1.229 nm, 
b = 0.305 nm, c = 0.581 nm, and j3 = 103.77° agree very well with the experimetally determined values and make this 
explanation rather unlikely. In order to get a deeper understanding of this energy shift, a precise description of the DF 
by means of model dielectric functions including an analysis in dependence on the temperature of the corresponding 
oscillators is necessary which is beyond the scope of this paper. Despite this energy shift between the theoretical and 
experimental DF we can, however, note that otherwise the structure of the DF and in particular also the anisotropies 
are reproduced perfectly. We observe the same relative shifts and ordering of the peaks and onsets for different light 
polarisations in Fig. 3a and also the peak heights are reproduced correctly. This means essentially, that the blue shift 
of the experimental spectra with respect to the theoretical ones is approximately a rigid shift. 

From the calculated electronic structure and the symmetry of the bands the polarization and magnitude of the 
transition from the highest valence bands into the lowest conduction bands at the T-point is derived. The involved 
valence bands which contribute mainly to the DF and the polarization of the corresponding transitions are also 
indicated in Fig. 3b. Except of the transitions from the first and second topmost valence bands into the conduction 
band which are dipole-allowed for a polarization mainly parallel to the z- and x-axis, respectively, the transitions are 
polarized either along the y-axis or within the x — z-plane. This is in agreement with the DF determined by means 
of ellipsometry where the observed transitions above the absorption edge in the tensor elements e xx and e zz occur at 
the same energy and differ from those which appear in £ yy . 


IV. SUMMARY 


To summarize, all tensor elements of the dielectric function of / 3 -Ga 203 were determined by spectroscopic ellipsom¬ 
etry and calculated by modern many-body theory including quasiparticle effects in the band structure calculations 
as well as excitonic and local-field effects in the calculation of the dielectric tensor. They have been shown to be 
in excellent agreement with each other. By means of the band structure calculation, the observed transitions in the 
experimental dielectric function were assigned to the corresponding valence bands which explains the observed energy 
difference in the onset of the absorption for a dipole polarized along the symmetry axis and within the x — 2 -plane 
of about 500 meV. Both, the experiment as well as the theoretical calculations, reveal a non-vanishing off-diagonal 
element which reaches values up to \e xz \ ~ 0.30 in the absorption range. This element is responsible for the rotation 
of the dielectric axis and although it is quite small, it causes an rotation of the dielectric axes of up to 20° in the 
transparent spectral range {\e xz \ < 0.02) . 
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